Monitoring Entanglement Evolution and Collective Quantum Dynamics 
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We generalize a recently developed scheme for monitoring coherent quantum dynamics with good 
time-resolution and low backaction [Reuther et at, Phys. Rev. Lett. 102, 033602 (2009)] to the case 
of more complex quantum dynamics of one or several qubits. The underlying idea is to measure with 
lock-in techniques the response of the quantum system to a high-frequency ac field. We demonstrate 
that this scheme also allows one to observe quantum dynamics with many frequency scales, such 
as that of a qubit undergoing Landau-Zener transitions. Moreover, we propose how to measure the 
entanglement between two qubits as well as the collective dynamics of qubit arrays. 
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I. INTRODUCTION 

A most fascinating question in quantum mechanics is 
how much a measurement acts back on a quantum sys- 
tem and thus influences the state of the latterji"— Ac- 
cording to the classic measurement postulate, the pro- 
jective measurement of a quantity causes a collapse of 
the wave function into an eigenstate of the correspond- 
ing operator. However, quantum mechanics only allows 
probabilistic statements about the state into which the 
wave function collapses. This evolution into a mixture 
can be modeled by coupling the quantum system via the 
measurement operator to a measurement apparatus, i.e., 
to a macroscopic environment which provokes dissipation 
and decoherence. A natural interpretation of this pro- 
cess is that information about the quantum system leaks 
into the environment, while the system acquires entropy. 
Since the environment is perceived as classical object, one 
may assume that the pointer of a measurement appara- 
tus corresponds to a collective environment coordinatci^ 
Then the question arises how the quantum system and 
this collective coordinate influence each other. In other 
words, which information about the quantum system is 
contained in the collective coordinate and how strong is 
the backaction of the environment on the system? 

A specific example for such a measurement device is 
a low-frequency "tank circuit" coupled to a supercon- 
ducting qubit i^i^ Here, one makes use of the fact that 
the resonance frequency of the oscillator depends on the 
qubit state which, in turn, influences the phase of the 
oscillator response. This allows one to measure both the 
charge and the flux degree of freedom of superconducting 
qubits. The drawback of this scheme, however, is that the 
coherent qubit dynamics is considerably faster than the 
driving. Thus, one can only observe the time-average of 
the qubit state, but not time-resolve its dynamics. Mea- 
suring the qubit directly by driving it at resonance is 
possible as wellii This however induces Rabi oscillations 
and, thus, alters the qubit dynamics significantlyi^i^. 

In Ref. lid , we proposed to probe a quantum system 



with a weak high-frequency drive acting directly upon 
the qubit. We found that the reflected signal possesses 
a time-dependent phase shift which is related to a qubit 
observable. This relation was validated for a Cooper- 
pair box that is in a superposition of the two lowest en- 
ergy eigenstates which form a qubit. In particular, it was 
shown that coherent oscillations of this qubit are visible 
in the reflected signal, while the driving-induced backac- 
tion stays at a tolerable level. This enables time-resolved 
monitoring of the coherent qubit dynamics. 

In this article we demonstrate that this measurement 
scheme is applicable also to more complex quantum dy- 
namics, and even the evolution of entanglement between 
two qubits can be observed in this way. In Sec. [Hi 
we introduce the underlying system-bath model and de- 
scribe how we quantify measurement fldelity and backac- 
tion. Moreover, we derive a relation between the reflected 
ac signal and the monitored time-dependent expectation 
value of the quantum system. This relation, which con- 
stitutes our measurement scheme, is numerically tested in 
Sec, mil for the case of a qubit undergoing Landau-Zener 
sweeps. In Sec. IIV[ we present a protocol for monitor- 
ing entanglement between two qubits. Finally, we derive 
in Sec. |V] how the collective dynamics of a qubit array 
enhances the signal. 



II. HIGH-FREQUENCY RESPONSE OF A 
QUANTUM SYSTEM 

A. System-bath model 

We consider a quantum system interacting with a dis- 
sipative environment, as depicted in Fig. [TJ The system- 
bath Hamiltonian is given byii"— 
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FIG. 1: (Color online) Sketch of the input-output formal- 
ism: A quantum system described by a Hamiltonian Hq is 
coupled via its charge degree of freedom Q to an environment 
in form of a transmission line with ohmic effective impedance 
Zq. When the system is probed by an input signal ^in, the 
output signal ^out contains information about the bare system 
dynamics. 



Here, Hq denotes the system Hamiltonian and Q is the 
system operator that couples to the environment. To be 
specific, we assume that Q is the system excess charge. 
In the realm of circuit QED, the bath is typically mod- 
eled by a transmission line in terms of LC circuits with 
charges qk and conjugate fluxes ipk, where Cfc and Lk 
are effective capacitances and inductances, respectively. 
Furthermore, Wfe = {LkCk)~^^^ denotes the angular fre- 
quency of mode qt, and Xk are the corresponding cou- 
pling constants in units of frequencies. The system-bath 
interaction is fully characterized by the spectral density 
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which is assumed to be ohmic with an effective impedance 
Zq, i.e. I{uj) = ujZqw^— Assuming a weak system-bath 
interaction, we obtain the Bloch-Redfield master equa- 
tion for the reduced system density operator p by stan- 
dard techniques,— ii^ 

P^CoP^ -'j\H,uP] - j:[Q, [Q,P]] - i^[Q, [Q.p]+], 

(3) 
where [A, B\+ = AB + BA denotes the anticommutator, 

Q = i[HQ,Q]/h and 



Q = - 

TT 



dr I duj S{uj)cos{ujt)Q{-t). (4) 

"'0 



Here, S{uj) = I{uj) coth{huj/2kBT) denotes the Fourier 
transformed of the symmetrically ordered equilibrium 
bath correlation function at temperature T. The short- 
hand notation X{t) represents the interaction-picture op- 
erator UQ{t)XUQ{t), where Uq denotes the system prop- 
agator. 



B. Input-output formalism 

The central concept of our measurement scheme is to 
relate the quantum dynamics of the externally probed 
central circuit to a response via the transmission line. 
In order to quantify the response, we employ the input- 
output formalism] ^^1^^ We start from the Heisenberg 



equation of motion (operator label H) for the environ- 
mental mode fc, 



gf + Lulq^ = XkiOkQ" , 



(5) 



which is linear in the bath operators. Its formal solution 
for initial time to = reads 



Ik W = QkiO) cos{u}kt) 




(/?fe(0)sin(wfet) 



+ Xk dTsm[cjk{t-T)]Q"{T). (6) 
Jo 

It depends on the history of the system, which gets ex- 
pressed through the convolution integral in the second 
line of Eq. (jH]). Inserting the obtained solution into the 
Heisenberg equation of motion for any system observable 
X yields 

k 



X / dTsm[uk{t-T)]Q"{T)--[Q'',X"] 




qk(Q) cos{uJkt) + ipk{0) sm{ujkt) 



(7) 



At this point we define the operator 

^k 



erw-E^^-(^ 



gfe(O) cos(wfet) + ipkiO) sm{u}kt) 



.(8) 

which allows a convenient notation. It is fully determined 
by the bath correlation function and, following the above 
equation of motion, Eq. ([7]), only depends on the envi- 
ronment operators at initial time. Therefore, it can be 
interpreted as a signal entering via the transmission line, 
i.e. describes the input noise acting upon the system. 

After a partial integration of the third term of Eq. (O , 
the counter-term proportional to [{Q^)'^,X^] is can- 
celed. Furthermore, for the ohmic spectral density 
I{u}) = ujZo, Eq. ([7]) becomes local in time, and one ar- 
rives at the quantum Langevin equation^Sr— 
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[H,",X"]-^[Q",X"]Q 



-^[Q^,X^]^rit)- (9) 

It contains the input field Chi™('^) ^-nd an ohmic dissipative 
term ex ZqQ^ . 

It is also possible to express the quantum Langevin 
equation in terms of the outgoing fiuctuations. This can 



be achieved by solving Eq. ([5]) for the bath modes starting 
at a later time ti > t. The solution is the time-reversed 
counterpart of Eq. (j6]) and reads 



Ik (0 = qk{ti)cos[uk{t-ti)] + i -—ipk{ti)sin[ujk{t-ti)] 

V ^k 
rti 
-Xk dTSm[LUk(t-T)]Q"{T) . (10) 



According to Eq. ([8]) , we define the sum over all modes k 
with respect to the first two terms of Eq. PH)) as the out- 
going noise Coiit(0- ^^ contrast to Chi™(*)j it is influenced 
by the system at earlier times t < ti and, thus, contains 
information about the system dynamics J^ Proceeding as 
above, we obtain the time-reversed Langevin equation 



X« 



[<,X^] 



+ l^[Q^^X"]Q"-'-[Q\x"]iZS)- (11) 

Its characteristic is a negative damping term, i.e., it is for- 
mally obtained from the quantum Langevin equation Q 
via the substitution Zq — >■ —Zq. This corresponds to 
backward propagation in time, as Eq. (|lip only accounts 
for the environment at later times ii > t. 

The difference between the two Langevin equations re- 
lates the input to the output fluctuations as^ 



ciw-erw = -2^oQ^ 



2iZo 



[H^.Q"]- (12) 



The last equality has been obtained from the Heisenberg 
equation of the charge operator, Q^ = [ilh)[HQ ,Q^], 
which is independent of the environment since Q 
commutes with the system-bath coupling Hamiltonian. 
Equation (|12p determines the influence of the system on 
its environment and, thus, represents the central relation 
of the input-output formalism. At this stage, it is worth 
emphasizing the environment possesses many degrees of 
freedom which all acquire information about the quan- 
tum system. Therefore, the collective bath coordinates 
^in/out (^) ^'''-' classical in the sense that their expectation 
values can be interpreted as the outcome of a single mea- 
surement i^ 



C. Response to high-frequency driving 

The central idea of our scheme is to excite the sys- 
tem by a classical ac drive and to measure the result- 
ing system response. Physically, the driving enters like 
the quantum noise via the transmission line. Therefore, 
it can be modelled as a coherently and highly excited 
bath mode in the classical limit weakly coupled to the 
system.— This means that the input fluctuations are 
augmented by a deterministic ac contribution. Thus, 
they become £,init) = ^h^i^) + Acos{p.t) and, accord- 
ingly, ^out(i) = Cout(^) + Acos{^t). The deterministic 



terms also affects the input-output relation (|T2l) which 
now reads 



Ut{t) = er(i) + Acos{m) ^[<, Q"] . (13) 



In the corresponding expectation value, where we return 
to the Schrodinger picture for convenience, the fluctua- 
tion ^'^™ vanishes, such that, 



iUtit)) ^Acos{nt) 



2iZn 



{[Ho,Q])t 



(14) 



This expresses the expectation value of the output field 
in terms of the input field and a system expectation 
value, (. . .)t = tr[/5(i) . . .]. The latter contains informa- 
tion about the system dynamics and is the central quan- 
tity of interest. 

The quantum Langevin equations © and (jlip have 
been convenient for deriving the relations between the in- 
put and the output fields, Eqs. ([T2|) and ([TS]). For the ac- 
tual computation of the system expectation value on the 
r.h.s. of Eq. ([T?)) . by contrast, such stochastic, operator- 
valued equations are less practical. Therefore, we derive 
in addition an equivalent quantum master equation which 
not only serves for numerical computations, but will also 
provide an analytical high-frequency approximation. 

We start by noticing that a classical ac drive Acos(rit) 
coupled to the system charge Q corresponds to a Hamil- 
tonian QAcos{nt). Thus, the system Hamiltonian Hq 
has to be replaced by 



H{t) = Ho + QAcos{Qt) 



(15) 



A key requirement of our measurement scheme is that 
the driving must not significantly alter the system dy- 
namics. Thus, we assume that the amplitude A is suffi- 
ciently small, which allows us to treat the driving per- 
turbativelyi^ In doing so, we use the ansatz p{t) = 
Po(t) + Pi(t), where po(t) describes the dynamics of the 
system without ac driving, and pi(t) denotes the correc- 
tion to the system state. With this ansatz, the master 
equation reads 



m^ [co+ci{t)][po{t)+p,it)] 



(16) 



where the perturbative driving manifests in the Liouvil- 
lian 



Ciit)p{t) ^ --A[Q, p{t)]cosnt , 
h 

To lowest order in A, the correction pi obeys 



(17) 



pi{t) ^ £oPi{t) - j:A[Q, poit)]cosint) . (18) 

This linear inhomogeneous equation of motion can be 
solved formally in terms of a convolution between the 
propagator of the undriven system and the inhomogene- 

ity, 



Pi{t) = e^«Vi(0) - ^y dTe^°(*-^U[Q,po(T)]cos(r!T) , 



(19) 



with pi(0) = 0. If the driving frequency fl is much larger 
than aU relevant system frequencies, the integral can be 
simplified by time-scale separation: First, we split the 
above integral into a sum of integrals over complete pe- 
riods and a final one over the remaining time. Assum- 
ing that the slow po{t) is practically constant during one 
oscillation period of 27r/f2, the integrals over complete 
driving periods vanish, while the last contribution can 
be evaluated using po{t) w po(t). In this way, we obtain 
the solution for the time evolution of the reduced density 
matrix 

p{t) = po{t) - -^ [Q, po{t)] sin{nt) , (20) 

where eA/hi^ is identified as the necessarily small per- 
turbation parameter. Within this approximation for the 
density operator, the expectation value of the output, 
Eq. P^ . becomes 



(CoutW) = Acos(m) 



2iZn 



[Ho,Q])o, 



+ -^{[[Ho,Q],Q]),^,sinm, (21) 

where the subscript in (. . .)o.t = tr[po(0 ■ • ■] refers to the 
undriven dynamics. 

At this point we note that (^out(i)) contains both low- 
frequency components stemming from the pure system 
dynamics as well as high-frequency components induced 
by the external driving. The latter corresponds to the 
second term of Eq. ()20[) . Now, in an experiment, it is 
feasible to single out the high-frequency components of 
the outgoing signal using a lock-in technique, where the 
incoming signal represents the reference oscillator. This 
removes the second term on the r.h.s. of Eq. ((2T|). such 
that the outgoing signal becomes 



(CtW>=^cos(m) 



2AZo 



[[Ha,Q],Q])^tsHm 



Since the second term is a perturbative correction, it has 
to be small so that Eq. ((22|) can be written as 



(23) 



{^lUt)) ^ Acos[m - cji,{t)] 
with the time-dependent phase shift 



2^0 



t>Ut) = h^mo,Q].Q] 



(24) 



This relation constitutes the basis of our measurement 
scheme. It connects a phase shift 0iJf(t) between the 
high-frequency input and the output signals to the low- 
frequency dynamics of a system observable. In other 
words, Eq. (|24p enables time-resolved monitoring of an 
open quantum system by measuring the phase shift (f)^^ [t) 
with lock-in techniques. We emphasize that our mea- 
surement scheme is rather generic and can in principle 
be applied to any open quantum system. 



It is important to note that there are situations in 
which the phase shift (i>]^f{t) is independent of the sys- 
tem dynamics or even vanishes. The latter is obviously 
the case when Q commutes with the system Hamilto- 
nian. A constant phase shift is obtained for an harmonic 
oscillator i/o = huj^la^ a + 1/2) that couples via its posi- 
tion, i.e., Q oi a + a\ as one easily obtains by evaluat- 
ing the double commutator in Eq. (|24p . This is evident 
from the superposition principle which holds for a linearly 
driven harmonic oscillator. However, already a coupling 
Q oc [a + a^Y is sufficient for obtaining a non-trivial 
response. 



D. Measurement fidelity and baclcaction 

As broached above, the relation between the phase of 
the output signal and a system observable relies on a 
high-frequency approximation. In an experiment, the 
driving frequency is finite, however. Thus, the ex- 
perimentally obtained phase (j\^^{t), which is actually 
recorded by the lock-in amplifier, may differ from the the- 
oretically predicted phase (/'hf(t) oc ([[ifo, Q]^Q])o,t- Dur- 
ing lock-in amplifying, the high-frequency components of 
(^out(O): ^'^d with it (?!)out (^)j ^-re extracted from the out- 
put signal (^out(O)j where the classical incoming signal 
(6n(0) = Acos{^t) serves as the reference oscillator. 

Thus, it is crucial to test with a numerical simulation 
how well both phases agree in a realistic case. As a crite- 
rion, we employ the measurement fidelity which we define 
as the normalized overlap 



F 



=(Cu?,0l!f) 



dt[Cu?W]' dt[<j^Ut)] 



-1/2 



dt 



C?(i)0Sf(O 



(25) 



with time integration over the duration of the measure- 
ment. The ideal value F = 1 corresponds to the case 
of perfect proportionality between the measured phase 
C?(0 and <f(i) a ([[iJo, Q], QDo,*- Note that cg,{t) is 
computed in the absence of the ac driving, while </>out (^) is 
obtained in its presence. Moreover, we have to mimic the 
action of the lock-in amplifier numerically, as described 
in Sec, mil in the second paragraph after Eq. (|5^ . 

Furthermore, for any quantum measurement, one has 
to worry about backaction on the system in terms of deco- 
herence. In our measurement scheme, decoherence plays 
a particular role, because both the driving and the ohmic 
environment couple to the central system via the same 
mechanism. This is reflected by the fact that the pre- 
dicted phase of Eq. (pi)) is proportional to the dissipa- 
tion strength Zq. In terms of the generalized, dimen- 
sionless dissipation strength a ~ e'^Zo/h, it is required 
that a < 0.1 in order to preserve a predominantly co- 
herent time evolution. In this context, our measurement 



is weak, but yet destructive in the sense that it rehes 
on the naturally provided interaction with a dissipative 
environment causing decoherence. This marks a signifi- 
cant difference to conventional quantum non-demolition 
(QND) measurement schemesiii^ Here by contrast, the 
additional decoherence due to the driving is minor. This 
is obvious from the fact that within first-order approxi- 
mation the system purity is not affected. 



Tr(p2) = Tr(pg) 



'0/' 



(26) 



which follows readily from Eq. ([20| . Only in this sense, 
our measurement scheme can be considered as being of 
non-demolition character. 

As a concrete measure for how much the driving per- 
turbs the system, we use in our numerical investigations 
the time average D of the trace distance^ 



1 



i?(f) = -Tr|p(<)-poWI 



(27) 



In the ideal case, D vanishes, while D = 1 if both density 
operators are completely unrelated. 



III. MONITORING LANDAU-ZENER SWEEPS 

In Ref. liO, we investigated the quality of our measure- 
ment scheme using the dynamics of a decaying qubit state 
as an elementary example. Here, we go a step further and 
demonstrate that it is applicable to more complex coher- 
ent dynamics. 

As an example, we consider a Landau-Zencr tunnel- 
ing process between the states of a charge qubit pre- 
pared in a Cooper- pair box (CPB)^i^ a paradigmatic 
example of complex quantum dynamics in a seemingly 
simple system. In particular, it can be used for state 
preparation2Si22 and entanglement generatio n^^i^^ in a 
qubit-resonator system. In recent experiments^!^ mul- 
tiple Landau-Zener sweeps were performed with a charge 
qubit coupled to a low-frequency tank oscillator whose 
resonance frequency depends on the time-averaged qubit 
state. This enabled the measurement of the accumulated 
phase which determines the qubit state. Here by con- 
trast, we demonstrate that a time- resolved observation 
of the complex dynamics during a single Landau-Zener 
transition is possible as well. 

We employ the CPB Hamiltonian 

H^^^ = 4Ec{N - N,r - - ^ (|iV+l)(7V|+h.c.), 

N=-oo 

(28) 
where N is the number of excess Cooper pairs in the box. 
The charge operator reads Q = 2eN = 2eY,N ^W)i^\ 
with e being the elementary charge. The charging en- 
ergy Ec is determined by the various capacitances of the 
CPB, while the scaled gate voltage Ng and the effective 
Josephson energy 6 are controllable. If the charging en- 
ergy is sufficiently large, and iV^ approaches a charge 



degeneracy point determined by the half-integer value 
^dog ^ j\j-^ _ \I2^ for some integer number TVo, only the 
two lowest charge states |iVo — 1) = \\.) and |A^o) = It) 
matter and form a qubit, while the energy gap to the 
higher charge states is much larger than the qubit split- 
ting.— li^i^ In the two-level approximation, the qubit is 
described by the Hamiltonian 



Hn ~ — ecTr (5(Jt 

2^2 



(29) 



The Pauli matrices Ui arc defined in the qubit subspace 
and e = 8Ec{No — Ng), where Ng ranges in the inter- 
val [Nq — l,-/Vo], is the effective qubit bias. Moreover, 
Qqb = ecTz while by virtue of relation (|24)) the phase of 
the output signal is linked to the qubit observable a^ 
according to 



,,. _ 4.e'ZoS 



o,t 



(30) 



Thus, the high-frequency component of (Q), which is 
manifest in the phase of the outgoing signal (fT2|) . contains 
information about the low-frequency qubit dynamics in 
terms of the unperturbed {ax)o- It is worth mentioning 
that the preceding discussion equally applies to the case 
of a flux or phase qubit with an inductive coupling to the 
transmission line. 

The qubit undergoes a Landau-Zener transition when 
the gate voltage in terms of the background charge Ng 
in the Hamiltonian (|28p is swept from a large positive 
to a negative finite value through the charge degeneracy 
point at N^°>^ with constant velocity, Ng{t) = N° + Vgt. 
Accordingly, the qubit Hamiltonian p9p depends on time 
and can be written as 



Hoit) 



v{t-to) 



-0-^ - -O-x 



(31) 



where we have defined v = SEcVg and to = ~N^/vg. 
Hamiltonian pip is a valid approximation on condition 
that the sweep ends timely prior to a subsequent anti- 
crossing in the CPB spectrum. 

Even if the qubit approximation describes the CPB 
faithfully for large charging energies, Ec ^ S, one must 
take into account excitations of higher charge states by 
the ingoing signal; see also the discussions in Ref. [lO- In 
our simulations, we operate at A^o = 6 and truncate the 
CPB Hilbert space at A^s = 12 states, which turned out 
to be sufficient to reach numerical convergence. 

Figure ^a,) illustrates the time evolution of the CPB 
observable 



Q 



ieS 



^ {\N){N + 1\-\N + 1){N\) 



(32) 



JV=-c 



during a Landau-Zener transition between the two lowest 
states; see inset of Fig. [D^a). Within two- level approxi- 
mation, this operator reads Q = {—eS/h)ay. The param- 
eters are chosen such that the final populations of the 
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FIG. 2: (Color online) Monitoring the dynamics of a dissipa- 
tive Landau-Zener sweep of a qubit, which is initiated in a 
weakly probed CPB with A^s = 12 states operated at A'^o — 6. 
The selected anticrossing lies at Ng = No — 1/2 = 5.5, and the 
initial qubit state is |iV = 5) = |t), while ||) = \N = 6). The 
other parameters are Ec = 2.6(5, Vg = 0.125 S/h, A = S/e, 
a = 0.05, and fi = 50 5/h. (a) Time evolution of the output 
signal {Q)t <x {£,out{t)) — {£,in{t)) of the full CPB, and its lock-in 
amplified phase (/f>out (i) (frequency window AQ = 165), com- 
pared to the estimated phase (^hf(i) oc {ax)o,t in the qubit 
approximation. Fast oscillations with frequency fl = 505/h 
stemming from the input signal are barely resolved. The in- 
set shows the three lowest CPB eigenenergies as a function 
of Ng{t). (b) Power spectrum of (Q) for the driven systems, 
comparing the full CPB Hamiltonian (solid) and in two-level 
approximation (dashed) . 



two states are roughly equal, because in this crossover 
regime between the adiabatic and the non-adiabatic dy- 
namics, coherent oscillations are most pronounced. Note 
that since we plot the time derivative of the charge. 



ay), the final population of the upper state. 



(Q) 

[1 -|- {(Tz){t — > oo)]/2, given by the Landau-Zener proba- 
bility Plz = exp(— 7r(5^/2?iu) sa 0.546, is not part of the 
information depicted in Fig. [2Ja). 

The backaction of the weak driving signal on the sys- 
tem dynamics in terms of a weak modulation of {Q)t 
is barely visible. The associated spectrum in terms of 
((5)tj, depicted in panel (b), nevertheless reflects the driv- 
ing in terms of a small peak at the driving frequency 
VL = 505/h. Due to the time-dependent qubit splitting, 
the spectrum exhibits a broad range of frequencies cen- 



tered at zero, reflected in the sideband structure around 
r2. In the time domain these sidebands correspond to a 
signal (^hf^(t)) = Acos[m - C?(0]- 

As is indicated above in Sees. Ill Cj and HTdI the phase 
^out (0 '^^ii be retrieved by lock- in amplification of the 
output signal in an experiment. We mimic this proce- 
dure numerically in the following wayi^ We only consider 
the spectrum of (^out(i)) in a frequency window 51 ± AfJ 
around the driving frequency and shift it by —17 in order 
to center it at zero. The inverse Fourier transformation 
of the spectrum cutout to the time domain yields the 
time-dependent phase 0out (0 which is expected to agree 
with (/)°f(i); see the discussion in Sec. Ill Dl According to 



Eq. (p4| . this phase reflects the unperturbed time evolu- 
tion of {(Tj:)o with respect to the qubit. 

Figure[2lja) shows that the numerically extracted phase 
'/'out (0 is indeed in a very good agreement with </'hf(0 oc 
{ax)o.t for an appropriate choice of parameters. This 
holds even if the condition of high-frequency probing, 
i.e. fl exceeding all relevant system frequencies by far, 
is not strictly fulfilled. Appropriate values for n need 
to be determined from the width of the sideband dis- 
tributions. These can be estimated from the duration 
tlz of a Landau-Zener transition which is tlz ~ 20h/S 
here. Consequently, the spectral peak has the width 
vTY,7,/h ~ 206 /h. Thus, in order to avoid any overlap 
between both distributions, the driving frequency must 
fulfill il > 206 /h. In the present case, Eqs. ^ and ^ 
are corroborated already for H, « 506 /h. 

Furthermore, we require small dissipation, a < 0.1, to 
keep the time evolution predominantly coherent. In this 
context, one needs to take care of the interplay of deco- 
herence and the external sweep velocity and its impact 
on the LZ probability and dynamics j^i^ii^^ 

Towards the end of the sweep, one notices a devia- 
tion between the predicted phase 0hf(i) and (/>out (0 which 
becomes curved upward. This stems from the presence 
of higher charge states, which get slightly populated as 
soon as Ng reaches the next half- integer value N^°^ + 1, 
that is, the subsequent avoided level anticrossing in the 
CPB spectrum. Accordingly, as visible in Fig. [DJb), 
the spectrum {Q)ui contains more pronounced frequency 
components than the one of the two-level approximation 
(Qqb)w However, the deviations between both spectra 
are minor, which corresponds to the good agreement be- 
tween 0ou?(t) and (plfit). 

The measurement fidelity defect 6F — 1 — F between 
(p'^^it) and 0hf(i), computed with Eq. (PSJ) . is depicted 
in Fig. inja) as a function of the driving frequency. Ap- 
parently, resonant excitations to non-qubit CPB states 
play a minor role here since, in contrast to the case of 
a time-independent qubit energyj^^ we do not observe 
characteristic resonance peaks in 5F. Since all CPB fre- 
quencies are varied in time during the LZ sweep, the pop- 
ulation of non-qubit states by the driving is mainly sup- 
pressed. Furthermore, the fidelity only slightly increases 
after the driving frequency n has exceeded a minimal 
value of vTi^z/h. This corroborates that $7 w 505/h is a 
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FIG. 3: (Color online) (a) Fidelity defect 5F = 1 - F and (b) 
time-averaged trace distance D between the density operators 
of a driven and an undriven CPB during a LZ sweep. Both 
quantities are depicted for various driving amplitudes A as a 
function of the driving frequency SI. All other parameters are 
the same as in Fig. (2] 



good choice. 

With regard to the driving amphtude, one has to find 
a compromise. The phase contrast of the outgoing sig- 
nal ([21]) obviously increases with A, whereas the driv- 
ing perturbs more and more the low-frequency dynamics. 
Here, the choice oi A = S/e yields a maximum fidelity of 
F ~ 0.9 for the above driving frequency. This apparently 
low value essentially stems from the discussed deviation 
between (j>'^{t) and 0[Jf(t). Far from a subsequent LZ 
sweep, the system is faithfully described with only the 
qubit levels. All in all, for these values, eA/hfl « 0.01- 
0.1, which justifies our perturbative treatment. 

The small-dissipation condition a ;< 0.1 together with 
the above conditions on the driving amplitude and fre- 
quency provides phase shifts of roughly (l\^ « 0.1°, 
which is small but still measurable with present tech- 
nologies. In particular, the constraint of large driving 
frequencies 57 for monitoring a LZ sweep reduces the max- 
imum visibility of the oscillations encoded in the phase 
shift. In contrast, a smaller value of J7, if applicable, may 
result in a visibility of the order of 1°. We turn back to 
these quantitative issues in Sec. IVl where we discuss the 
collective dynamics of equal quantum systems. 

Generally, decoherence may be influenced by ac driv- 
ing i^ Here, however, we do not observe a significant 
change of decoherence. We conclude this from Fig. ^h) 
which depicts the time-averaged trace distance, com- 



puted with Eq. (j27|) as a function of the driving frequency 
for various driving amplitudes. For intermediate ampli- 
tudes, we find roughly D oc A/fl, unless A is very small, 
which agrees with relation (P5| . This states that at high 
frequencies, the driving does not add decoherence, which 
confirms the picture drawn by investigating the measure- 
ment fidelity F. 

Our model Hamiltonian does not consider the excita- 
tion of quasi-particles in the superconductor. Therefore 
it is valid only as long as the driving frequency stays be- 
low the gap energy. Thus, a setup made of aluminum 
the driving is limited to fl < 100 GHz. With a typical 
Josephson energy of the order of some GHz, the driv- 
ing frequency required to monitor the dynamics of a LZ 
sweep already comes close to this limit. An implemen- 
tation with niobium, whose energy gap is considerably 
larger— and which has been used in circuit QED experi- 
mentsi^l should be less critical. 



IV. TWO-QUBIT ENTANGLEMENT 

When considering a second qubit, it may be entangled 
with the first one. In the present context, this raises the 
question whether our measurement scheme is sensitive 
to this property. In order to dwell into this question, we 
focus on a system of two charge qubits described by the 
Hamiltonian 



i?2c0 = !^(ai^) 



^i^O+S'xi^V, 



(1)^(2) 



(33) 



where the upper indices label the qubits. Here, we as- 
sume that both qubits possess the same energy splitting 
Lu and are mutually coupled with strength g. A qubit 
coupling scheme as described here is effectively realized 
in the case of, e.g., two CPBs within two- level approxi- 
mation which are capacitively coupled to a transmission 
line resonator. In analogy to the above discussions, our 
measurement scheme can be applied on condition that 
the driving field is not in resonance with transitions to 
higher energy levels of the underlying physical system J^ 
If the qubit-resonator interaction is dispersive, i.e., if the 
resonator is far detuned with respect to the qubits, it 
mediates an effective qubit-qubit XX-couplingr^ such 
that the two-qubit system can be described by the effec- 
tive Hamiltonian (|33| . Alternatively, an effective XX- 
coupling between two qubits can be mediated by a third, 
dispersively detuned qubitj^ 

Both charge qubits couple capacitively to a common 
environment via their charge operator 

g2q = e(aW+af). (34) 

An adequate measurement of the qubit-qubit entangle- 
ment is the concurrence 

C = max{xi - X2 - X3 - X4, 0} , (35) 

where the Xj denote the ordered square roots of the eigen- 
values of the matrix p{(Jy ay )p*{cry (Jy ) with p denot- 
ing the two-qubit density matrix4£ 



For later convenience we express the system operators 
in terms of the four, maximaUy entangled Bell states, 



|a>±) = -i=(|oo)±|ii)) 

|*±) = -i=(|10)±|01)) 



(36) 

(37) 



In the basis {|$+), |$_), 1*+), |vI/_)} the two-qubit 
Hamiltonian (1331) reads 



■ff2q,o - ;: 



/ .g -w \ 

-uj -g 

5 

V -gj 



(38) 



while the charge operator becomes 
Q2q = 2e 



1 0^ 

10 



oy 



(39) 



This representation of the coupling to the environment 
evidences that neither the system-bath coupling nor the 
ac driving generate transitions between the subspaces 
spanned by the states {|$±)} and {|*±)}, respectively. 
Concerning the latter states, it is even such that they 
are not coupled to the environment at all, i.e., they 
form a decoherence-free subspace4ii^ For this reason, 
we henceforth focus on initial preparations in the sub- 
space spanned by {|<i>±)}. In this particular case, the 
concurrence is given by 



C, 



{|*± 



.}(0 = 2|pooii(t)| = 2|Tr{p(t)|00)(ll|}|, (40) 



i.e., by an off-diagonal density matrix element in the 
product basis. In terms of the pseudo-spin operators a'^ 
and a'^ defined in the basis {|$+), |<I>_)}, the Hamilto- 
nian and the charge operator become 



Qaq = 2ecr^ 



(41) 
(42) 



respectively. 

At this point, we recall Eq. (P^ which relates the phase 
of the outgoing signal to a particular observable of the un- 
driven system. Inserting the subspace Hamiltonian and 
charge operator, Eqs. ([4T|) and (|42p. we obtain 



(t) = ^ ia. 



Ml 



z/0,t ' 



(43) 



with the dimensionless damping strength a = c'^Zq/H, 
and 

{<j',)o,t ^ Tr{po(t)^;} - 2Rc{p°„n(i)} • (44) 

In other words, the phase of the outgoing signal is pro- 
portional to the real part of the density matrix element 
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FIG. 4: (Color online) Monitoring the entanglement evo- 
lution of two coupled qubits with initial state 1$+). The 
lock-in amplified phase cj)^ (blue dashed line) is compared 
to the estimated phase (f'hii^) (r^d sofid line) which is pro- 
portional to the real part of the density matrix element 
Pooii(^)- The latter is needed for reconstructing the concur- 
rence C{t) = 2|pooii(f)| of the unperturbed system. Parame- 
ters are Q, — Ibcj, g = 0.1 uj, A — 0.1 oj and a — 0.08. 



Pooii(0 of the unperturbed system. Figure U demon- 
strates that this relation holds sufficiently well: It com- 
pares the numerically calculated phase of the lock- in mea- 
surement, 0Q^j (i), to the phase (j>1[(t) ex (o'^)o,t predicted 
by the measurement relation 1^^. We find an excellent 
agreement between both quantities, while the maximum 
angular visibility reaches values of around 1°. Thus, we 
find that it is possible to access directly the key density 
matrix element needed to reconstruct the concurrence 
C{|$^)}(t) for the subspace {|$±)}. 

Monitoring R.e{pooii(*)} already provides a lower 
bound for the entanglement and may even provide a good 
estimate for it. Nevertheless, it is desirable to measure 
also its imaginary part. For this purpose one has to ro- 
tate the basis of the second qubit by applying the phase 
gate 



u: 



(2) 



phase 



1 
i 



(45) 



This local unitary transform does not alter the two- 
qubit entanglement at the time it is applied. For small 
values of 5, it approximately commutes with the to- 
tal time evolution operator of the system, cxp{iHt/h). 

[exp{iHt/h),U^y^'^^J = for g = 0. In analogy to 
(UH); the relation for the phase shift then becomes 



I.e. 
Eq, 



0hf(t)«MPooiiW} 



(46) 



Thus, we have access to both the real and the imaginary 
part of the density matrix element Pooii(i)- One just 
needs to initiate two copies of a pure entangled state in 
the subspace {|$±)}, apply the gate of Eq. P5)) to one of 
them, and monitor their dissipative time evolution. The 
difference between the two-qubit dynamics in either case 



(2) 

due to the phase gate C/phasc i^ relevant only at large times 
when the qubits approach their stationary state. For 
the transient stage depicted in Fig. HI we have verified 
that the phase of the output signal indeed reflects the 
imaginary part of the concurrence (not shown). This, 
in turn, enables one to reconstruct the characteristics 
of the concurrence C{|$^)i.(i) and thus of the amount 
of entanglement between the two qubits, according to 
Eq. (|40l). 

In a similar way, one finds access to the entanglement 
when the system is initially prepared in the subspace 
\'^±). This can be done by applying the so-called Pauli-X 
gate to the second qubit. 



a 



(2) 



r(2) 



1 

1 



(47) 



effectively leaving us in the |<&±) subspace again, from 
where one can pursue as described above. Like this, it is 
possible to monitor the time characteristics of entangle- 
ment of all kinds of entangled two-qubit states, as long 
as the subspaces {|<&±)} and {|*±)} are not mixed. 

Admittedly, the relation between the entanglement 
and 0°f, manifest in Eqs. (P5)) and (pi)) . holds only if the 
concurrence can be expressed as the expectation value of 
an observable providing the required off-diagonal matrix 
element of the density operator. Obviously, this docs not 
always apply. However, two qubits prepared in a Bell 
state represent such a case which is of relevance and in- 
terest. 

So far, we have addressed initially pure entangled 
states, which constitutes an idealization and is not al- 
ways achievable in an experiment. A more realistic sce- 
nario copes with the drawback that one starts from a 
mixed state which naturally possesses less entanglement. 
Moreover, the concurrence is then no longer determined 
by a single off-diagonal matrix element [cf. Eq. (|40p ]. but 
rather has to be computed from the definition of the two- 
qubit concurrence, Eq. pSI) . In order to test how much 
information about the concurrence C(t) is nevertheless 
contained in the measurement signal (j331)> we study the 
time evolution for the preparation in an unpolarized mix- 
ture or Werner state^^ 

pw = {l-p)\<i>+){<^+\ 

+ I (|$_)($_| + |*+)(*+| + |*-)(*-|) , (48) 

where p denotes the degree of depolarization. A perfect 
initial preparation in the fully entangled state l^^-) is 
then characterized by p = 0. For p > 0.5, the state is 
separable. In the realm of circuit QED, a polarization 
degree oi p = 0.1 has already been achieved;^ while in 
a quantum optical context, even the much lower value 
p < 0.01 could be reached4^ 

Our goal is now to compare the concurrence C(t) with 
C/|$_|_\}(i) which holds for preparation in a Bell state. 
As a measure for the agreement, we use again the fidelity 
F defined as the scaled overlap between both quantities. 
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FIG. 5: (Color online) Fidelity defect SF = 1 - F of the 
concurrence curves C'{|#±)}(i) and C{t) in dependence of the 
degree of depolarization, p. The function overlap is integrated 
over the decay time of the two-qubit state. As discussed in 
the text, the upper bound (5Knax = 0.03 marks the largest 
acceptable value of p. All other parameters as in Fig. [l] 



i.e., by Eq. (PS|) but with the phases replaced by the con- 
currences. Figure [5] shows the resulting fidelity defect 
6F = 1 — F as a function of the initial degree of depolar- 
ization. From a numerical comparison of the concurrence 
curves C{t) and C'{|$j^)}(i) for different values oi SF (not 
shown), we set the threshold for the maximum accept- 
able fidelity defect to 5i^max = 0.03. The investigation of 
several examples showed that, beyond a fidelity defect of 
this order, both curves no longer appear similar. Thus, it 
appears that the measured phase shift captures the con- 
currence dynamics satisfactorily only for a SF below this 
threshold. Remarkably, the latter is not reached until an 
initial depolarization of p < 0.1. From this, we conclude 
that our measurement protocol allows one to gain infor- 
mation about qubit-qubit entanglement in circuit QED 
even for depolarized initial states. 



V. COLLECTIVE QUANTUM DYNAMICS 

The typical signal that can be reached with our proto- 
col, as regards measurements of single quantum systems, 
is a phase shift of the order 1°, as estimated in Sec. IIIII 
Of practical interest is whether this already marks the 
maximum reachable. In the following, we suggest how 
to increase the visibility by monitoring the collective dy- 
namics of a large number of identical quantum objects 
at the same time. For this purpose, we exemplarily in- 
vestigate an "atomic cloud" containing a large number 
of atoms as described in Ref. |4y. There, an ensemble of 
cold atoms interacting with a superconducting cavity, de- 
scribed by the Dicke model, was investigated. Provided 
that it is possible to single out one particular atomic tran- 
sition, the cloud represents an ensemble of identical two 
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level atoms. This gives rise to the A'^-qubit Hamiltonian 



VI. CONCLUSIONS 



N 



Hm 






N 



(49) 



The interaction between the electromagnetic field in the 
cavity and the atoms is assumed dipole-like and can be 
cast into the operator QNqia + o,^)- Here, the collective 
coupling operator QNq is chosen dimensionless and given 



N 



QNq = / ,0- 



U) 



NS. 



(50) 



and Sz = J2i=i '^z' /\N denotes the collective spin of 
the cloud. 

In the case where an open transmission line is em- 
ployed in the setup rather than the (closed) cavity, a 
continuum of modes replaces the single cavity mode in 
the modeling. In this sense, the transmission line can 
be understood as one-dimensional heat bath. The dipole 
interaction between the cloud and the environment be- 
comes Hint = QNqJ2k Cfe(al + flfc)- Thus, one naturally 
ends up with the Caldeira-Legget model of Eq. ^ where, 
formally, Ck = XkUJk/Ck- 

Assuming an ohmic environment, the phase shift (j24p 
then becomes 



€fit) 



2a 



[[HoMM)^, 



(51) 



where a denotes the dimensionless damping strength, 
and a5 is the coupling strength to the environment per 
atom. Inserting Hpfq.Q and QNq into the measurement 
relation ([5T|) yields 



Kfit) = 



iaSVN 



\^x)o,i 



(52) 



with Sx ~ X]i=i cTx /ViV. Thus, the phase signal scales 
with ^/N as compared to Eq. ([29| . By virtue of this 
relation, monitoring the collective dynamics of TV two- 
level systems via the total spin paves the way towards a 
higher angular visibility. 

Since the validity of the measurement relation (|5ip is 
restricted to small angles, the visibility may be enhanced 
up to the order 10°. In order to be specific, we evaluate 
(j)^f for the parameters'^ d/2TTh = 6.8GHz, a = 6 • 10-^ 
ri/27r = 20GHz. This yields an effective single-atom cou- 
pling to the ohmic environment of a6/2Trh sa 0.4MIIz. 
Despite this very small system-bath interaction strength 
per atom, the estimated phase resolution for N = 10^ 
atoms in the cloud reaches (p^f = 4° — 5°. At the same 
time, atomic decoherence is drastically reduced. Under 
these conditions, which are experimentally achievable,— 
the real-time monitoring of an atomic ensemble consti- 
tutes a major improvement in angular visibility in com- 
parison with the single-qubit case. 



We have derived a scheme for monitoring the low- 
frequency dynamics of a quantum system, where a partic- 
ular focus has been put on quantum circuits. The central 
idea is to couple the quantum system capacitively to a 
high-frequency ac field and to measure the reflected sig- 
nal. We have shown that the phase of the latter contains 
information about a particular expectation value of the 
central system. In an experiment, this information can 
be extracted with lock-in techniques. This scheme has 
some similarities to the qubit readout via a tank circuili^'^ 
or a quantum point contact 4i The essential difference is 
that the proposed high-frequency driving enables a time- 
resolved measurement. 

The underlying relation (j24p between a particular ex- 
pectation value and the phase shift of the reflected signal 
relies on a time-scale separation. This constitutes an ap- 
proximation, which has to be validated. In so doing, 
we have demonstrated that the measurement fidelity is 
rather good. As an example, we have studied a qubit un- 
dergoing Landau-Zener transitions. Owing to the time 
dependent level splitting of the qubit, its dynamics pos- 
sesses a broad frequency spectrum. Thus, it represents 
a more challenging test case than the coherent oscilla- 
tions considered in Ref . Hol. Also for the present case, the 
measurement fidelity turned out to be rather good. 

A further relevant issue is the backaction of the mea- 
surement on the system. The relevant backaction is de- 
coherence which here is unavoidable because the driving 
and the environmental degrees of freedom couple to the 
system via the same mechanism. This is reflected by the 
fact that the phase shift being the recorded signal is pro- 
portional to the decoherence rate. Thus a stronger cou- 
pling to the driving can be achieved only at the expense 
of more decoherence. 

Concerning possible applications in solid-state quan- 
tum information processing, the generalization towards 
systems with two or more qubits is rather important. For 
the case of two qubits, we have demonstrated that it is 
possible to monitor the time evolution of their entangle- 
ment, provided that one can access the relevant density 
matrix elements via the phase shift of the reflected sig- 
nal. Even though the entanglement decays due to the 
unavoidable decoherence, it is still possible to monitor 
an appreciable number of cycles between entangled and 
separable states of interacting qubits. Furthermore, our 
protocol produces reliable results even in the case where 
the monitored state is mixed rather than pure. 

Finally, we have also addressed the question of how 
much stronger the measurement signal can be if the sin- 
gle qubit is replaced by a whole array of qubits that un- 
dergo the same dynamics. The result is that the signal 
scales with the square root of the number of qubits. This 
even leaves some room for improving the fidelity and re- 
ducing the coupling to each qubit. As a consequence, 
decoherence becomes less relevant, which is encouraging 
in particular when attempting a proof-of-principle real- 
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ization. 



Acknowledgments 

We thank Jens Siewert and Enrique Solano for fruit- 
ful discussions. This work has been supported by 



DFG through the Collaborative Research Center SFB 
631, project A5, and through the German Excellence 
Initiative via "Nanosystems Initiative Munich (NIM)". 
G.R. acknowledges financial support from the NIM seed 
funding program. D.Z. acknowledges financial support 
from MICINN under grant numbers FIS2008-01240 and 
FIS2009-13364-C02-0. 



Electronic address: 'georg.reuther@piiysik.uni-augsburg.de 
^ D. F. Walls and G. J. Milburn, Quantum Optics (Springer, 

Heidelberg, 1995), 2nd ed. 
^ M. A. Nielsen and I. L. Chuang, Quantum Computing and 

Quantum Information (Cambridge University Press, Cam- 
bridge, 2000). 
^ H. P. Breuer and F. Petruccione, The theory of open quan- 
tum systems (Oxford University Press, Oxford, 2002). 
* W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 
^ M. Grajcar, A. Izmalkov, E. Il'ichev, T. Wagner, 

N. Ouklianski, U. Hiibner, T. May I. Zhilyaev, H. E. 

Hoenig, Ya. S. Greenberg, et al., Phys. Rev. B 69, 

060501 (R) (2004). 
^ M. A. Sillanpaa, T. Lehtinen, A. Paila, Y. Makhlin, 

L. Roschier, and P. J. Hakonen, Phys. Rev. Lett. 95, 

206806 (2005). 
'^ Ya. S. Greenberg, E. Il'ichev, and A. Izmalkov, EPL 72, 

880 (2005). 
« Y.-X. Liu, C. P. Sun, and F. Nori, Phys. Rev. A 74, 052321 

(2006). 
^ M. Wilson, Phys. Today 60, 20 (2007). 
^° G. M. Reuther, D. Zueco, P. Hanggi, and S. Kohler, Phys. 

Rev. Lett. 102, 033602 (2009). 
" A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. 

Fisher, A. Garg, and W. Zwerger, Rev. Mod. Phys. 59, 

1 (1987). 
^^ P. Hanggi, P. Talkner, and M. Borkovec, Rev. Mod. Phys. 

62, 251 (1990). 
^■^ Y. MakhUn and A. D. Mirlin, Phys. Rev. Lett. 87, 276803 

(2001). 
" B. Yurke and J. S. Denker, Phys. Rev. A 29, 1419 (1984). 
^^ M. H. Devoret, Quantum Fluctuations in Electrical cir- 
cuits (Elsevier, Amsterdam, 1995), chap. 10, Les Houches, 

Session LXIII. 
^^ K. M. Fonseca-Romero, S. Kohler, and P. Hanggi, Phys. 

Rev. Lett. 95, 140502 (2005). 
" M. Grifoni and P. Hanggi, Phys. Rep. 304, 229 (1998). 
^® C. W. Gardiner and M. J. CoUett, Phys. Rev. A 31, 3761 

(1985). 
^^ G. Johansson, L. Tornberg, and C. M. Wilson, Phys. Rev. 

B 74, 100504(R) (2006). 
20 G. W. Ford, J. T. Lewis, and R. F. O'ConneU, J. Stat. 

Phys. 53, 439 (1988). 
2^ P. Hanggi, F. Marchesoni, and F. Nori, Ann. Phys. 

(Leipzig) 14, 51 (2005). 
22 P. Hanggi and F. Marchesoni, Rev. Mod. Phys. 81, 387 

(2009). 
2-' S. Kohler, T. Dittrich, and P. Hanggi, Phys. Rev. E 55, 

300 (1997). 
2* G. Wendin and V. S. Shumeiko, in Handbook of Theoretical 



and Computational Nanotechnology, edited by M. Rieth 

and W. Schommers (American Scientific Publishers, Los 

Angeles, 2006), vol. 3, p. 223. 
2^ S. N. Shevchenko, S. Ashab, and F. Nori, Phys. Rep. 492, 

1 (2010). 
2'' K. Saito, M. Wubs, S. Kohler, P. Hanggi, and 

Y. Kayanuma, EPL 76, 22 (2006). 
2'^ D. Zueco, G. M. Reuther, P. Hanggi, and S. Kohler, Phys- 

ica E 42, 363 (2010). 
2* M. Wubs, S. Kohler, and P. Hanggi, Physica E 40, 187 

(2007). 
2^^ D. Zueco, P. Hanggi, and S. Kohler, New J. Phys. 10, 

115012 (2008). 
^° W. D. Oliver, Y. Yu, J. C. Lee, K. K. Berggren, L. S. 

Levitov, and T. P. Orlando, Science 310, 1653 (2005). 
^^ M. Sillanpaa, T. Lehtinen, A. Paila, Y. Makhlin, and 

P. Hakonen, Phys. Rev. Lett. 96, 187002 (2006). 
^2 D. Vion, A. Aassime, A. Cottet, P. Joyez, H. Pothier, 

C. Urbina, D. Esteve, and M. H. Devoret, Science 296, 

886 (2002). 
^^ J. H. Scofield, Am. J. Phys. 62, 129 (1994). 
^'' P. Nalbach and M. Thorwart, Phys. Rev. Lett. 103, 220401 

(2009). 
^^ M. Wubs, K. Saito, S. Kohler, P. Hanggi, and 

Y. Kayanuma, Phys. Rev. Lett. 97, 200404 (2006). 
^^ V. Novotny and P. P. M. Meincke, J. Low. Temp. Phys. 

18, 147 (1975). 
^■^ D. M. Berns, M. S. Rudner, S. O. Valenzuela, K. K. 

Berggren, W. D. Oliver, L. S. Levitov, and T. P. Orlando, 

Nature 455, 51 (2008). 
^* D. Zueco, G. M. Reuther, S. Kohler, and P. Hanggi, Phys. 

Rev. A 80, 033846 (2009). 
^® A. O. Niskanen, K. Harrabi, F. Yoshihara, Y. Nakamura, 

S. Lloyd, and J. S. Tsai, Science 316, 723 (2007). 
■*° W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
■*^ D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. 

Lett. 81, 2594 (1998). 
''2 T. Yu and J. H. Eberly Phys. Rev. B 66, 193306 (2002). 
"^ R. F. Werner, Phys. Rev. A 40, 4277 (1989). 
"'"' M. Ansmann, H. Wang, R. C. Bialczak, M. Hofheinz, 

E. Lucero, M. Neeley, A. D. O'ConneU, D. Sank, M. Wei- 

des, J. Wenner, et al.. Nature (London) 461, 504 (2009). 
""^ P. Lougovski, E. Solano, and H. Wahher, Phys. Rev. A 71, 

013811 (2005). 
'^^ J. Verdii, H. Zoubi, C. KoUer, J. Majer, H. Ritsch, and 

J. Schmiedmayer, Phys. Rev. Lett. 103, 043603 (2009). 
^"^ S. Ashhab, J. Q. You, and F. Nori, New J. Phys. 11, 083017 

(2009). 



